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Zusammenfassung 

Neues Anwachsen der italienischen Wolfspopulation: eine genetische Untersuchung mittels 
Mikrosatelliten. 

Um die genetische Differenzierung innerhalb der anwachsenden italienischen Wolfspopulation ( Ca¬ 
m's lupus) zu untersuchen, wurden 38 Individuen von 4 verschiedenen Stellen des derzeitigen Ver- 
breitungsgebietes mittels Mikrosatelliten analysiert. Die 6 betrachteten Mikrosatelliten loci zeigten 
einen hohen Grad an Polymorphismus und die kombinierte Identitatswahrscheinlichkeit reichte 
von 1(T 4 bis 10" 6 . Insgesamt war die gefundene DNA-Variation betrachtlich und kaum niedriger 
als jene beim nordamerikanischen Grauwolf. Die beiden altesten italienischen Teilpopulationen, 
jene in der Toskana und jene im zentralen bis sudlichen Apennin, zeigten nur eine maftige Diver- 
genz. Die gefundenen Werte stimmten mit einer Abstammung der ersteren von der Westalpen-Sub- 
population uberein, wahrend die letztere Ahnlichkeit zur Subpopulation der Westkiiste zeigte. Der 
Genfluft innerhalb des italienischen Wolfsbestandes war hoch. Die geographische Verteilung der 
Probengebiete stimmte mit den genetischen Abstanden iiberein, was auf das Vorliegen von isola¬ 
tion durch Distanz' hindeutet. Hohe Divergenzniveaus wurden zwischen den italienischen Proben 
und jenen aus anderen europaischen Gebieten gefunden. Insgesamt deuten die Ergebnisse darauf 
hin, daft, trotz der Einformigkeit im Bereich der mtDNA, der italienische Wolf noch immer einen ho¬ 
hen Grad an Kern-DNA-Variabilitat und eine ausgepragte genetische Identitat besitzt. Eine weitere 
Vergrofterung des Wolfsareals, wie sie aufgrund der hohen Dispersionsrate zu erwarten ist, konnte 
in einigen Jahrzehnten dazu fiihren, daft die historische Isolation von einem genetischen Aus- 
tausch mit dem osteuropaischen Genpool abgelost wird. 
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Abstract 

Climatic fluctuation led to isolation of some populations of temperate species on southern moun- 
taintops during warming trends. The most southern species of arvicoline rodents ( Microtus guate- 
malensis, M. oaxacensis, M. quasiater, and M. umbrosus) in North America may be relicts, isolated 
in the mountains of Mexico and Guatemala at the end of the Pleistocene. We used parsimony and 
likelihood analyses of complete mitochondrial cytochrome b gene sequences of 28 species of Micro¬ 
tus, including eight Eurasian species, holarctic M. oeconomus, and all extant North American spe¬ 
cies except the island endemic M. breweri. North American species of Microtus were monophyletic 
under the maximum-likelihood criterion, but polyphyletic under parsimony. Likelihood ratio tests 
and bootstrapping indicated a rapid basal radiation with short intervals between cladogenic events. 
However, several sister taxon relationships were robust to bootstrapping or consistent between 
methods. We found that M. quasiater was sister to M. pinetorum, and these taxa were sister to a 
clade of M. oaxacensis and M. guatemalensis. The phylogenetic position of M. umbrosus, however, 
was unclear. Monophyly of the relicts was rejected by a likelihood ratio test, suggesting multiple 
southern invasions by arvicoline rodents. Phylogenetic data for these and other co-distributed taxa 
should be used in conservation efforts for these remote areas. 

Key words: Microtus, Arvicolinae, historical biogeography, conservation, molecular systematics, 
Mexico 


Introduction 

Climatic fluctuations have been linked to 
the spatial expansion and contraction of 
species (Graham and Grimm 1990; Graham 
et al. 1996). Isolation during contraction 
phases may have stimulated allopatric di¬ 
versification, particularly at the southern 
edge of ranges. For example, the mountains 
of Mexico and Guatemala host a highly en¬ 


demic flora and fauna (Ramamoorthy et 
al. 1993) that includes a diverse set of or¬ 
ganisms associated with mesic environ¬ 
ments. These organisms apparently invaded 
the region during cooler periods and then 
became isolated at higher elevations as con¬ 
ditions at lower elevations became drier 
and warmer. This invasion and isolation 
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cycle may have occurred numerous times 
during the Pleistocene and has contributed 
to the complex biogeographic history of 
the region (Sullivan et al. 1997). 

The genus Microtus is holarctic and north 
temperate in distribution and in the New 
World reaches its southern limit in Central 
America. At higher latitudes, up to five spe¬ 
cies (e. g. M. longicaudus, M. miurus, M. 
oeconomus, M. pennsylvanicus , and M. 
xanthognathus in Yukon Territory) may be 
found in close proximity. However, species 
tend to be more allopatrically distributed at 
southern latitudes. For example, M. guate¬ 
malensis, M. oaxacensis, M. quasiater, and 
M. umbrosus are endemic to separate 
mountains in the cloud and pine forests of 
Mexico and Guatemala (Fig. 1). M. quasia¬ 
ter , the Jalapan vole, is found in the southern 
Sierra Madre Oriental in Central Mexico. 
M. oaxacensis , the Tarabundi vole, is iso¬ 
lated in the Sierra de Juarez of Oaxaca, and 
occurs at elevations of 1,600 (Sanchez et 
al. 1996) to 2,499 m (Jones and Genoways 
1967). M. umbrosus , the Zempoaltepec vole, 
is restricted to approximately 80 km 2 at ele¬ 
vations ranging from 1,829 to 3,000 m (Frey 
and Cervantes 1997) around Mt. Zempoal¬ 
tepec in Oaxaca. M. guatemalensis , the most 
southern species, occurs from the mountains 
of central Chiapas south to central Guate¬ 
mala. These species may be the result of 
peripheral isolation of ancestors that were 
more widely distributed during cool periods 


of the early to middle Pleistocene (Hoff¬ 
mann and Koeppl 1985). 

In contrast to these four disjunct species, M. 
mexicanus is widespread in Mexico and oc¬ 
curs in limited sympatry or parapatry with 
M. oaxacensis, M. quasiater , and M. umbro¬ 
sus with a discontinuous distribution ex¬ 
tending from New Mexico and Arizona to 
southern Mexico (Hall 1981). Its fossil re¬ 
cord is limited to the late Wisconsinan from 
San Josecito in northeastern Mexico and lo¬ 
calities farther north (Zakrzewski 1985). 
Phylogenetic relationships between M. 
mexicanus and other species of Microtus 
are unclear, but there is little indication that 
this species shares a common ancestor with 
other Mesoamerican species. Thus, its pre¬ 
sence in Mexico probably reflects an inde¬ 
pendent colonization. 

The aim of our study therefore is to investi¬ 
gate the phylogenetic relationships of the 
southern species of Microtus in order to test 
hypotheses regarding their ancestry in 
North America and their taxonomic rela¬ 
tionships to other species. 


Material and methods 

DNA was extracted from ethanol-preserved tis¬ 
sues of these four southern Microtus (Tab. 1) with 
methods described in Conroy and Cook (1999). 
The cytochrome b gene (hereafter cyt b) was am¬ 
plified in two or three sections and sequenced in 


Table 1. Specimens of Meso-American species of Microtus examined in this study (in addition to those reported 
in Conroy and Cook [2000]) including location from which each specimen was collected, and CNMA (Coleccion 
Nacionalde Mamiferos, Universidad Nacional Autonoma de Mexico) catalog number. 


Species 

Collection location 

CNMA Catalog Number 

M. umbrosus 

Mexico: Oaxaca Cerro Zempoaltepetl, 5 km N Sta Ma. Yacochi. 
Mpio. Tlahuitoltepec, 2450 m. 

34890, 34894 

M. guatemalensis 

Mexico: Chiapas: Cerro Tzontehuitz, 13 km NE San Cristobal 
de las Casas, Mpio. Chamula, 2880 m. 

35262 


Mexico: Oaxaca: 11 km SW La Esperanza, Mpio. Santiago 


M. oaxacensis 

Comaltepec, 2000 m.; Oaxaca: 11 km SE La Esperanza, Mpio. 
Santiago Comaltepec, 2000 m. 

27415, 33815 

M. quasiater 

Mexico: Veracruz: 5 km W Naolinco, Mpio. Naolinco, 1650 m. 

35282,35274 
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both directions (Perkin-Elmer Prism® dye termi¬ 
nator kit; Fst-RR, 402119) on an ABI 373a auto¬ 
mated sequencer. We included cyt b sequences 


from 24 species of Microtus and two species of 
Clethrionomys (Conroy and Cook 2000). In the 
phylogenetic analysis we represented each species 



Fig. 1. Map of the distribution of species of Microtus in North America and south into Guatemala (redrawn from 
Hoffmann and Koeppl 1985). 
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with a single individual, but we examined intra¬ 
specific variation where possible by including 
multiple representatives for 24 of the 28 species. 
We used unweighted parsimony (MP) and maxi¬ 
mum likelihood (ML) with the software PAUP*, 
version 4.0d64 (Swofford 1998). We estimated 
parameters for likelihood models of increasing 
complexity (JC: Jukes and Cantor 1969; 
HKY85: Hasegawa et al. 1985; and GTR: Yang 
1994 a; HKY85 + L, and GTR + T: Yang 1994 b). 
We tested these with likelihood ratio tests to de¬ 
termine significant differences among their likeli¬ 
hood scores. We used 1,000 random addition se¬ 
quences to locate multiple tree islands in the 
heuristic MP searches. As tests of node strength, 
we bootstrapped the parsimony analysis 1,000 
times and bootstrapped the final ML analysis 100 
times. We tested several alternative phylogenetic 
topologies against the ML tree by forcing likeli¬ 
hood searches to find the ML topology with parti¬ 
cular constraints and then comparing scores 
(Kishino and Hasegawa 1989). We tested mono- 
phyly of the four relictual Mesoamerican species, 
and the four shortest maximum parsimony trees. 
Besides examining relationships among these 
southern latitude taxa, several tests were con¬ 
ducted to ascertain the effects of adding these 
four species to more general hypotheses concern¬ 
ing the systematics of Microtus previously con¬ 
ducted (Conroy and Cook 2000). We tested the 
monophyly of subgenus Stenocranius (Rausch 
1964). This subgenus includes the Asiatic species 
M. gregalis and the North American M. miurns 
and M. abbreviatus. It has been hypothesized that 
the species in this subgenus diverged relatively re¬ 
cently. Secondly, we tested the monophyly of a 
clade of taiga adapted voles ( M . richardsoni, M. 
xanthognathus, M. chrotorrhinus) that was not 
supported in our previous study (Conroy and 
Cook 2000). And, third, we tested the monophyly 
of all North American taiga voles (Hoffmann 
and Koeppl 1985), because previously we were 
able to reject this hypothesis with the Kishino- 
Hasegawa test. By including these tests, our inves¬ 
tigation explores whether the addition of four 
southern taxa changes our interpretation of arvi- 
coline history in North America. Prior to this 
study, the relationships of southern species to 
other North American species was unclear. 

Results 

Base composition and distribution of vari¬ 
able sites of the cyt b gene was similar to 
other Microtus (Conroy and Cook 2000) as 


well as other mammals (Irwin et al. 1991). 
Of the 1143 base pairs, 472 were variable 
across 28 species of Microtus and two species 
of Clethrionomys. When outgroups were ex¬ 
cluded, 460 sites were variable. Of these, 
100 were in the first position, 23 in the sec¬ 
ond position, and 337 in the third position 
of codons. Of the 381 amino acids, 76 (20%) 
were variable across species of Microtus and 
the replacement pattern was consistent with 
structural models (e. g. Irwin et al. 1991). 
There were 351 parsimony informative nu¬ 
cleotide sites and gi statistics (gi = -0.325) 
indicated phylogenetic signal in the data. 
Maximum parsimony searches recovered 
four equally parsimonious trees (Fig. 2), 
each including a basal clade of the North 
American M. ochrogaster and Asian M. gre¬ 
galis. A clade of M. oeconomus, M. midden- 
dorffi, M. montebelli, M. kikuchii, and M. 
fords (hereafter the ‘Asian clade”), and 
the M. pennsylvanicus clade (i. e. M. penn- 
sylvanicus , M. montanus , M. townsendii , 
and M. canicaudus) were present in the 
four trees. M. pinetorum and M. quasiater 
were sister in all trees, had high bootstrap 
support in MP and ML analyses (99%) and 
relatively high decay values (12). M. oaxa- 
censis and M. guatemalensis were sister taxa 
in three of four MP trees, but the branch 
leading to this pair had weak bootstrap sup¬ 
port (<50%). Other clades were found in 
three or four of these shortest trees, but 
bootstrap support was generally low across 
basal relationships. M. umbrosus stemmed 
from the polytomy at the base of the Euro¬ 
pean and North American species. Random 
addition of taxa uncovered these four short¬ 
est trees, but simple addition (i. e. alphabe¬ 
tical by species name) led to an island of 
longer length trees. 

In the ML analysis, the HKY85 + T likeli¬ 
hood model (transition/transversion ratio = 
3.43, a - 0.213) was chosen since more com¬ 
plex models (e. g. GTR + T) were not 
significantly more likely and produced the 
same topology (not shown). As in other stu¬ 
dies of cytochrome b (e. g. Sullivan et al. 
1997), the addition of the gamma-distribu¬ 
ted rate parameter, allowing among-site 
rate variation, contributed significantly to 
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Fig. 2. Strict consensus of four trees from the MP search. Each tree had a length of 2038, a consistency index of 
0.331, and retention index of 0.360. Values to the left of slash are bootstraps percentages greater than 50% from 
1,000 searches with simple addition of taxa; values to the right of slash indicate decay indices calculated with 
10 random-addition replicates for each search. 
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the likelihood. This model produded one 
tree (Fig. 3) in which M. gregalis was basal, 
followed by the Asian clade. North Ameri¬ 
can endemic species formed a clade and 
the European species formed a sister clade 
to North American species. The Mesoamer- 
ican endemics were not basal within the 
clade of North American species. Three 


Mesoamerican species (Microtus guatema- 
lensis, M. oaxacensis, and M. quasiater) dis¬ 
played the same branching pattern as three 
of the parsimony trees, while M. umbrosus 
was sister to M. chrotorrhinus. Other rela¬ 
tionships were similar to previous analyses 
(Conroy and Cook 2000). For example, the 
M. Pennsylvaniais and Asian clades were 
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Fig. 3. Maximum likelihood phylogenetic tree based on the HKY85 + T model (see text for parameter values). 
Values above branches are bootstrap percentages greater than 50% from 100 bootstraps. 
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supported and M. mexicanus was sister to 
M. californicus. In Kishino-Hasegawa tests 
(Tab. 2), few of the alternate topologies 
could be rejected. Those that were rejected 
were one of the MP trees and two ML trees 
constrained to monophyly of the four 
southern species. 

Relative depth of divergence was estimated 
with pairwise likelihood distances between 
taxa estimated under the same model used 
for the ML phylogeny. M. abbreviatus and 
M. miurus , which differ by only 0.015, prob¬ 
ably separated at the end of the Pleistocene 
when rising sea levels isolated M. abbrevia¬ 
tus on islands in the Bering Strait. Another 
late Pleistocene divergence may be M. cani- 
caudus and M. townsendii which differ by 
0.058. Other pairs that were relatively clo¬ 
sely related included M. arvalis and M. ros- 
siaemeridionalis (0.072) and M. montanus 
and M. pennsylvanicus (0.086). Divergence 
values between the southern species (M 


quasiater and M. pinetorum [0.094], M. gua - 
temalensis and M. oaxacensis [0.113], and 
M. chrotorrhinus and M. umbrosus [0.137]) 
are greater than any of the preceding exam¬ 
ples, suggesting older speciation events. 
These comparisons should be interpreted 
with caution. Demographic and historical 
differences among species can alter rates of 
evolution and species of Milrotus probably 
have evolved at different rates. Also, genet¬ 
ic differentiation may not necessarily coin¬ 
cide with vicariant events, such as the rise 
of the Bering Sea. Finally, molecular esti¬ 
mates of divergence, particularly for a sin¬ 
gle locus, have substantial amounts of sto¬ 
chastic variation (Ayala 1999). 

Discussion 

Previous investigations into the history of 
North American species of Microtus sug- 


Table 2. Kishino-Hasegawa tests of tree topologies. The "Stenocranius” constraint forced M. miurus, M. gregatis, 
and M. abbreviatus to be monophylyetic. The "Second taiga vole" constraint enforced M. xanthognathus, M. chro- 
tonrhinus, and M. nchardsoni monophyly. The "All taiga voles" constraint enforced the monophyly of the "Second 
taiga vole" species with a clade of M. pennsylvanicus, M . montanus, M. townsendii, and M. canicaudus . The 
"Meso-Americans" constraint enforced M. guatemalensis, M. oaxacensis, M. quasiater, and M. umbrosus mono¬ 
phyly. Two trees were obtained from the ML search under this last constraint. The parsimony trees were derived 
from unweighted heuristic searches with ten random-addition replicates. Columns below indicate, from left to 
right, 1) -log likelihood score of that tree, 2) difference in score between that tree and the best tree, 3) the stan¬ 
dard deviation of that comparison, 4) the associated T test statistic, and 4) the associated p value. 


Tree 

Constraint 

-log likelihood 

Diff-lnL 

s. d. (diff) 

T 

P* 

1 

None 

9868.811 

(best) 




2 

Stenocranius Monophyletic 

9891.258 

22.447 

16.723 

1.342 

0.180 

3 

Second taiga vole clade 

9879.176 

10.365 

9.080 

1.142 

0.254 

4 

All Taiga voles Monophyletic 

9897.044 

28.233 

18.904 

1.494 

0.136 

5 

Meso-Americans monophyletic I 

9910.346 

41.535 

11.613 

3.577 

0.0004* 

6 

Meso-Americans monophyletic II 

9910.346 

41.535 

11.613 

3.577 

0.0004* 

7 

Parsimony Tree 1 

9908.109 

39.298 

22.586 

1.740 

0.082 

8 

ParsimonyTree 2 

9908.743 

39.932 

20.492 

1.949 

0.052 

9 

ParsimonyTree 3 

9908.493 

39.682 

20.581 

1.928 

0.054 

10 

ParsimonyTree 4 

9919.839 

51.027 

20.528 

2.486 

0.013* 


Significant at P < 0.05 
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gested an evolutionary history closely tied 
to fluctuating boreal ecosystems (Hoff¬ 
mann and Koeppl 1985). The elevational 
distribution of species of Microtus in Mexi¬ 
co and Guatemala indicates Holocene iso¬ 
lation due to the contraction of cool, moist 
forests. The paucity of synapomorphic mor¬ 
phologic characters, but apparent abun¬ 
dance of autapomorphic characters has led 
to their characterization as ancient, highly 
divergent species; however, their biogeo¬ 
graphic history has been enigmatic. This 
study reconsiders the biogeographic and 
evolutionary history of the southern relicts 
in light of molecular characters examined 
within a wider taxonomic sampling (28 spe¬ 
cies) for the genus. 

Biogeography 

The Mesoamerican species of Microtus were 
not basal in the North American clade, nor 
were they monophyletic. M. quasiater is clo¬ 
sely related to M. pinetorum , and M. oaxa- 
censis and M. guatemalensis are sister to this 
pair, while the fourth (M. umbrosus) is a 
weakly-supported sister taxon to M. chrotor- 
rhinus. Our data do not support the hypoth¬ 
esis of Hoffmann and Koeppl (1985) that 
these southern species of Microtus may be 
relicts of an early colonization of North 
America prior to the arrival of the ancestors 
of the other species that are endemic to the 
higher latitudes of North America. Interspe¬ 
cific distances do, however, suggest they be¬ 
gan to diverge earlier than suspected late- 
Pleistocene peripheral isolates such as M . 
abbreviatus and M. miurus . Assuming an 
equivalent rate of cyt b evolution across spe¬ 
cies of Microtus and Clethrionomys , Mesoa¬ 
merican Microtus apparently fragmented 
into isolated populations prior to the coloni¬ 
zation of North America by Clethrionomys , 
which is thought to have occurred during 
the early Pleistocene (Repenning et al. 
1990). The split between North American 
Clethrionomys gapperi and Eurasian C. glar- 
eolus (Tamura-Nei D = 0.079) occurred 
about the time their common ancestor in¬ 
vaded North America. However, disparate 
levels of ancestral polymorphism and his¬ 


tories of population fluctuation in these in¬ 
dependent colonizers, among other vari¬ 
ables, would impact estimates of ages of 
differentiation. Two topologies that con¬ 
strained ML searches for monophyly of 
these southern species were significantly less 
likely than the best ML topology, refuting 
common ancestry of Mesoamerican species 
of Microtus. 

Sullivan et al. (1997) and Sullivan et al. 
(2000) found that the Isthmus of Tehuante¬ 
pec was a significant geographic barrier for 
other mesic rodents such as Reithrodon- 
tomys and the Peromyscus aztecus group in 
Mexico and Guatemala. They summarized 
divergence patterns in other co-distributed 
organisms (see Ramamoorthy et al. 1993) 
and recommended that the significance of 
the isthmus as a biogeographic barrier be 
tested with other mesic taxa. M. guatema¬ 
lensis occurs east of the isthmus and is sister 
to M. oaxacensis , found west of the isthmus, 
suggesting that their common ancestor may 
have been distributed across the isthmus 
when conditions were cooler. Western and 
eastern populations subsequently diverged. 
Microtus, Reithrodontomys , and Peromys¬ 
cus are widely sympatric in this region, and 
levels of genetic differentiation across the 
isthmus are surprisingly similar with 
D = 0.075-0.091 in Peromyscus (Sullivan 
et al. 1997), D = 0.063-0.085 in Reithrodon¬ 
tomys (Sullivan et al. 2000), and 
D = 0.100 between M. guatemalensis and 
M. oaxacensis. If we assume similar rates 
of molecular evolution, these three clades 
may have diverged across the Isthmus of 
Tehuantepec at roughly the same time. 
High levels of biotic diversity in the region 
provide ample material to further test hy¬ 
potheses regarding the timing and number 
of mammalian colonizations of the region. 
The biogeography of other organisms, such 
as reptiles (e. g. Campbell 1984; Campbell 
et al. 1989), may also suggest former con¬ 
nections across the Isthmus of Tehuantepec. 

System a tics 

Microtus quasiater is a member of the sub¬ 
genus Pitymys (Musser and Carleton 
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1993) and shares dental morphology with 
extinct Microtus (Pitymys) meadensis , a 
widespread species of mid-Pleistocene 
North America and Mexico (Repenning 
1983). M. pinetorum has often been classi¬ 
fied as a member of Pitymys based on 
shared dental morphology (Zakrzewski 
1985). Thus, it is not surprising that mor¬ 
phological characters (Musser and Carle- 
ton 1993) and these DNA sequences place 
M. quasiater sister to M. pinetorum. M. qua- 
siater previously was considered sister to M. 
ochrogaster (Moore and Janecek 1990) in 
an allozyme study, but only nine of 19 North 
American species and no Palearctic species 
were examined. 

The evolutionary relationships of M. guate- 
malensis and M. oaxacensis have not been 
addressed in detail, although they have 
not been suspected to be closely related 
(Musser and Carleton 1993). Microtus 
guatemalensis is in the monotypic subgenus 
Herpetomys , but may have affinities with 
Pitymys (Martin 1987). Our analyses con¬ 
tradict its placement in Herpetomys , or 
suggestions that it is related to Phena- 
comys (Hinton 1926). The relationship be¬ 
tween M. oaxacensis and other species has 
also been obscure (Musser and Carleton 
1993), but it also was considered a part of 
an early pitymyine invasion (Hoffmann 
and Koeppl 1985; Martin 1974). A wide¬ 
spread ancestor (e.g. M. meadensis) may 
have given rise to M. pinetorum , M. guate¬ 
malensis , M. quasiater , and M. oaxacensis , 
prior to peripheral isolation in the eastern 
deciduous forests and southern cloud for¬ 
ests (Hoffmann and Koeppl 1985). The 
mitochondrial phylogeny suggests isolation 
occurred first between an ancestor of 
M. pinetorurn-M. quasiater and an ances¬ 
tor of M. guatemalensis-M. oaxacensis 
although the branching order among these 
clades is weakly supported. The latter pair 
may have diverged after invasion across 
the Isthmus of Tehuantepec, while the for¬ 
mer pair diverged following an episode of 
range retraction induced by climatic warm¬ 
ing. Significant portions of the temperate 
flora on mountain tops in Mesoamerica 
are as old as the late Cenozoic (Graham 


1999) , however, we do not know of a 
comparable phylogenetic assessment of 
vegetation. 

A sister relationship between M. umbrosus 
and M. chrotorrhinus was weakly supported 
in the ML tree and was unexpected because 
they differ morphologically. Microtus um¬ 
brosus is the sole member of the subgenus 
Orthriomys (Musser and Carleton 1993) 
and has been considered a relict from an 
early invasion from Asia during the mid- 
Pleistocene by the extinct Phaiomys (Mar¬ 
tin 1987). Though previously considered 
closely related to M. xanthognathus (Hall 
and Kelson 1959), M. chrotorrhinus was la¬ 
ter distinguished based on chromosomal 
complement (Rausch and Rausch 1974). 
The lack of similarity between M. chrotor¬ 
rhinus and M. umbrosus and the “pity- 
myine” species suggests an independent in¬ 
vasion of the southern latitudes by a 
common ancestor of M. umbrosus and M. 
chrotorrhinus. 

North American monophyly and other 
phylogenetic hypotheses 

Our assessment of Microtus includes 28 of 
approximately 65 species of Microtus (Mus¬ 
ser and Carleton 1993), with the addition 
of Mesoamerican taxa potentially contri¬ 
buting to more accurate estimations of rela¬ 
tionships (Hillis 1996). This larger study 
did not significantly alter the ML topology 
previously obtained based on 24 species of 
Microtus (Conroy and Cook 2000). Mono¬ 
phyly of all North American species in¬ 
cluded in this study and of the M. pennsyl- 
vanicus clades was supported. The 
expanded analysis suggested two differ¬ 
ences based on Kishino-Hasegawa tests. 
Without Mesoamerican species, we rejected 
the topology (p = 0.026; Conroy and Cook 

2000) which constrained all North Ameri¬ 
can taiga voles as monophyletic (sensu 
Hoffmann and Koeppl 1985). However, 
with the inclusion of Mesoamerican species, 
this hypothesis was not rejected (p = 0.136). 
The expanded analysis also rejected a to¬ 
pology placing M. chrotorrhinus basal to 
all species of Microtus sampled and M. gre- 


Mexican Microtus 


341 


galis (Russia) within a clade of North 
American species (p = 0.013). 
Morphological material for interspecific 
comparison is abundant, but the phyloge¬ 
netic utility of some morphological charac¬ 
ters for arvicolines (e. g. tooth pattern) has 
been criticized due to high variation within 
and between species (Guthrie 1965; Zakr- 
zewski 1985). Despite the availability of 
standard karyotypes for many species of 
Microtus , non-differentially stained chro¬ 
mosomes have minimal phylogenetic infor¬ 
mation because the rate of chromosomal 
evolution varies greatly among species 
(Cervantes et al. 1997; Modi 1987). Indeed, 
our phylogeny suggests a complicated series 
of events are needed to explain chromoso¬ 
mal rearrangements in Microtus. Species 
with low diploid numbers are not sister to 
each other. For example, the mitochondrial 
data indicate M. oaxacensis (2N = 30) is sis¬ 
ter to M. guatemalensis (2N = 52), and M. 
canicaudus (2N = 24) is sister to M. town- 
sendii (2N = 50). These relationships sup¬ 
port Cervantes et al.’s (1997) contention 
that the interpretation of chromosomal evo¬ 
lution based only on standard karyotypes 
will be difficult. 

Conservation 

Mexico has one of the richest mammalian 
faunas due partially to an overlap of neo¬ 
tropical and nearctic biomes (Fa and Mo¬ 
rales 1993). Conservation efforts are com¬ 
plicated by this high diversity and a variety 
of threats (Ceballos and Navarro 1991). 
Species of Microtus in the region are nearc¬ 
tic relicts that are dependent on high eleva¬ 
tion, mesic habitats. Protection of habitat in 
the mountains of Oaxaca, a region of high 
mammalian diversity (Arita et al. 1997), 
could help conserve three endemics, M. 
oaxacensis , M. umbrosus , and M. qnasiater. 
Whether conservation criteria focus on rar¬ 
ity, diversity, or degree of endemism, these 
southern relicts warrant conservation 
concern. 

Molecular systematic studies of other ende¬ 
mic taxa should be considered in planning 
conservation efforts in this region (Baker 


et al. 1995; Cook et al. 2001). Our analysis 
suggests that temporal scale may be a cru¬ 
cial component to interpreting the signifi¬ 
cance of biogeographic barriers. Detection 
of patterns at different temporal scales 
could help resolve shared histories of taxa 
in the region (Avise 1994) and be used to 
conserve historical associations of flora and 
fauna. 

Species of Microtus in Mexico and Guate¬ 
mala are not monophyletic but instead are 
the result of at least three colonizations of 
this region during the Pleistocene: one by 
the ancestors of M. oaxacensis , M. quasia- 
ter , and M. guatemalensis , a second by the 
ancestor of M. umbrosus , and a third colo¬ 
nization apparently gave rise to M. mexica- 
nus. A lack of fossils inhibits the dating of 
cladogenic events among these species; 
however, depth of divergence relative to 
other splits within Microtus suggest mid- 
Pleistocene divergence. Morphological si¬ 
milarity between extant species (e.g. M. 
pinetorum and M. quasiater) and formerly 
widespread (mid to late Pleistocene) taxa 
that are now extinct (e.g. M. meadensis) 
suggest isolation by range retraction is a 
viable hypothesis. Morphological studies 
also support the shared history of several 
of the pitymyine species. The weakly-sup- 
ported sister relationship between M. um¬ 
brosus and M. chrotorrhinus has not been 
predicted previously and should be tested 
further. Phylogenetic analysis of other Me- 
soamerican organisms may help identify 
regions of shared evolutionary history and 
the role of significant biogeographic bar¬ 
riers in promoting diversification in this 
biologically rich region. 
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Zusammenfassung 

Die phylogenetische Stellung sudlicher Reliktarten von Microtus (Muridae Rodentia) 
in Nordamerika. 

KLimatische FLuktuationen fuhrten in Zeiten der Erwarmung zur Isolation einiger Populationen von Ar- 
ten aus gemaGigten Zonen auf siidlichen Bergen. Die sudlichsten Arten arvicoliner Nager ( Microtus 
guatemalensis, M. oaxacensis, M. quasiater und M. umbrosus) in Nordamerika mogen Relikte darstel- 
len, die am Ende des PLeistozans in den Bergen von Mexico und Guatemala isoliert wurden. Wir unter- 
suchten mittels Parsimony- und Likelihood-Analysen die kompletten mitochondrialen Cytochrom 
b-Gen-Sequenzen von 28 Arten der Gattung Microtus. Darunter befanden sich acht eurasische Arten, 
holarktische M. oeconomus und alle heute lebenden Arten Nordamerikas aufter der inselendemischen 
M. breweri. Die nordamerikanischen Arten waren in der Maximum-Likelihood-Analyse monophyletisch, 
in der Maximum-Parsimony-Analyse jedoch polyphyletisch. Likelihood-Ratio-Tests und Bootstrap- 
Analysen wiesen auf eine rasche adaptive Radiation mit kurzen Intervallen zwischen Kladogenese- 
Ereignissen hin. Die Analysen von Schwestertaxa zeigten jedoch Konstanz in Bootstrap-Analysen oder 
unterschiedlichen phylogenetischen Auswertemethoden. Nach unseren Ergebnissen war M. quasiater 
die Schwestergruppe von M. pinetorum und beide bildeten eine Schwestergruppe zu einer Klade aus 
M. oaxacensis und M. guatemalensis. Die phylogenetische Stellung von M. umbrosus blieb unklar. Die 
Monophylie der Relikte erwies sich nach Likelihood-Ratio-Tests als unhaltbar, was auf multiple sud- 
liche Invasionen arvicoliner Nager hinweist. Die phylogenetischen Daten fur 17 diese und andere Taxa 
sollten in Arbeiten zur Arterhaltung in diesen abgelegenen Gebieten berucksichtigt werden. 
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